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Pj abstract 



We present ATCA results of a 3 and 7 mm continuum survey of 20 T Tauri 
stars in the Chamaeleon and Lupus star forming regions. This survey aims to identify 
protoplanetary discs with signs of grain growth. We detected 90% of the sources at 
3 and 7 mm, and determined the spectral slopes, dust opacity indices and dust disc 
masses. We also present temporal monitoring results of a small sub-set of sources at 
7, 15 mm and 3+6 cm to investigate grain growth to cm sizes and constrain emission 
mechanisms in these sources. Additionally, we investigated the potential correlation 
between grain growth signatures in the infrared (10 /im silicate feature) and millimetre 
(1-3 mm spectral slope, a). 

Eleven sources at 3 and 7 mm have dominant thermal dust emission up to 7 mm, 
with 7 of these having a 1-3 mm dust opacity index less than unity, suggesting grain 
growth up to at least mm sizes. The Chamaeleon sources observed at 15 mm and 
beyond show the presence of excess emission from an ionised wind and/or chromo- 
spheric emission. Long-timescale monitoring at 7 mm indicated that cm-sized pebbles 
are present in at least four sources. Short-timescale monitoring at 15 mm suggests 
the excess emission is from thermal free-free emission. Finally, a weak correlation was 
found between the strength of the 10 /^m feature and a, suggesting simultaneous dust 
evolution of the inner and outer parts of the disc. This survey shows that grain growth 
up to cm-sized pebbles and the presence of excess emission at 15 mm and beyond are 
common in these systems, and that temporal monitoring is required to disentangle 
these emission mechanisms. 

Key words: protoplanetary disc, stars: variables: T Tauri, radiation mechanisms: 
thermal, radio continuum: planetary systems 



1 INTRODUCTION 

Observational signatures of grain growth in protoplanetary 
discs can be obtained from both the infrared (IR) and mil- 
limetre (mm) bands, which probe different regions of the 
disc and are sensitive to different grain sizes. IR emission 
comes from the warm inner (~l-5 AU) and upper layers 
of the disc, while the mm emission comes from the cooler 
outer disc regions (> 10 AU) and mid-plane where the bulk 
of the dust resides - see Fig. [T] The smallest submicron- 
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sized grains are strongly coupled to the gas, but as they 
grow they begin to decouple from the g as and settle to the 
mid-plane, resulting in grain s ize sorting (|Chiane: et al .1200 ll : 
iDuUemond fc Dominik |2004 ). Since thermal dust emission 
is dominated by dus t grains whos e size is similar to the ob- 
serving wavelength (|Drainell2006f ). multi- wavelength obser- 
vations are therefore needed to understand the distribution 
of the dust phase in protoplanetary discs. 

The strength and shape of the 10 li ra silicate feature 
can b e used as a grain growth indicator. iPrzygodda et all 
(|2003D found that low mass T Tauri stars with submicron- 
sized grains typically have a strong triangular shaped 
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Figure 1. Schematic of a p rotop lanetary disc, adapted from 
iDullemond fc Monnieij l|201of l and iMuIders fc DominikI l|2012f l. 
The grey area shows the gas within one (dark) and two (Ught) 
pressure scale heights. The black areas show the approximate 
region of the infrared emission. The millimetre and centimetre 
emission, tracing the larger dust grains, originates mainly from 
the mid-plane of the outer disc. Note that the origin of the ob- 
served infrared emission strongly depends on the orientation of 
the disc. 



10 /im feature, while discs with micron-sized grains typi- 
cally have weaker a nd broader 10 ^,m featur es. Similar re- 
sults were found by Ivan Boekel et all (j2003h for their sur- 
vey of intermediate mass Herbig Ae/Be stars. An equiv- 
alent effec t was suggested t o be caused by crystallisa - 



. iHonda et all l20Q3l : iKessler-Silacci et all l2006h . 
iKessler-Silacci et al.l (j2Q06l ) found an increase in flux near 



tion (e.g., 



11.3 /xm can be caused by polycyclic aromatic hydrocarbon 
and/or crystalline forsterite, which can mislead the grain 
growth interpretation. 

The 1-3 mm spectral slope, ai_3mm (hereafter a), where 
the flux F oc i/" and v is the frequency, can be used to esti- 
mate the dust opacity index /3 where dust opacity n„ . 
Assuming the emission is optically thin, the dust opacity in- 
dex can be written as /3 ~ a - 2 and an opacity index /3 1 
indicates grain growth up to mm sizes jPrain 3 l2006l ). This 
approach has been used to detect mm-sized grains in ~50 



discs to date (jNatta et al 


|2004|; lAndrews & Williamsll2005l: 


iLommen et al. 20071. 200£ 


, 2010l: Ricci et al.ll2010al bl). How- 



ever, care must be taken when relating a and P, as shallow 
spectral slopes can result from an ext ended optically thin 
disc or a compact optically thick disc (jBeckwith fc Sargent! 
Il99l] ). To break this degeneracy the di scs have to be spat ially 
resolved at mm wavelengths. Recentlv lRicci et al.l (|2012h de- 
termined that although only a small fraction of optically 
thick material is required to create a shallow spectral slope 
in the disc where the maximum grain size « 0.1 mm, the re- 
quired dust overdensities are much larger than the physical 
processes which concentrate dust (e.g., streaming instabili- 
ties, gravitational instabilities) can achieve. This strength- 
ens the notion that shallow spectral slopes at long wave- 
lengths are due to large mm/cm-sized grains in the outer 
parts of the disc. 

By extending the dust opacity index relationship to 
7 mm and beyond, it is possible to determine if grains up to 



cm sizes are present in discs (iTesti et al.ll2003l: IWilner et al.l 
I2OO5I : iRodmann et al.ll2006l : iLommen et al.ll2009l '). However 
for wavelengths of 7 mm and longer, the emission can result 
from a range of physical processes in young stellar objects, 
including thermal emission from dust, free-free emission 
from an i onised wind, and non-thermal chromospheric emis - 
sion (e.g.. lDullemond et al.|[2007l : iMillan-Gabet et al.ll2007l '). 
While chromospheric emission will be unresolved without 
the use of very long baseline interferometry, both thermal 
dust emission from the disc and free-free emission from a 
wind may be extended, and so simply resolving the emission 
is not enough to determine the dominant emission process. 

The spectral slope from 7 mm and beyond can be used 
to break this degeneracy. Depending on optical depth, two 
spectral slopes are known for free-free emission from spher- 
ically symmetric, constant velocity, ionise d wind. For an 
opaq ue wind the spectral slope, as = 0.6 JPanagia fc Felli 



1975 ) , while for optically thin wind, as = -0.1 (iMezger et al. 
1967 ). Using these definitions. iRodmann et al.l (|2006l ) usee 



the cm spectral slope from 2 to 3.6 cm of four T Tauri stars 
to determine that ~20% of the 7 mm flux was due to free-free 
emission. 

Another way to determine the emission mechanisms 
present in a disc is by temporal monitoring of flux variability 
at 7 and 15 mm. Thermal dust emission is expected to be 
constant over time, while thermal free-free emission from an 
ionised wind is expected to vary by ~20-40% on a tim escale 
of years (|Gonzalez fc Cantd|[20o3 : iLoinard et al.ll2007l ). Non- 
thermal emission varies over a timescale of minutes to hours 
by an order of m agnitude or more (jKutner et al.l Il986l : 
IChiang et anil996D . 

The results from IR and mm surveys have been used 
to study a potential correlation between the 10 /jm silicate 
feature and the millimetre s p ectral slope (ILommen et al.l 
l2007l . I2OIOI : iRicci et al.ll2010bt ). ILommen et al.l (|2007l . I2OIOD 
suggested a tentative correlation between the strength and 
shape of the 10 ^m silicate feature and a for a sample of 
Chamaeleon, Lupus, and Taurus sources exists. A correla- 
tion between these grain growth signatures would suggest 
simultaneous dust evolution of the inner and outer parts 
of the disc. The submicron-sized grains would no longer 
be replenished as grains grow to mm to cm sizes, which 
flattens the 10 /im silicate feature, while the simultaneous 
growth of larger particles in the midpl anc leads to a shal- 
lower mrn_s£ectral_slo£e (|Lommen et al...2007. . ,2010i ). How- 
ever, iRicci et al.l (|2010lj ) found no such correlation in their 
sample of p Ophiuchus and Taurus sources, and suggested 
that the small sample size and some inconsistency between 
the spectral slopes of specific sources could have lead to the 
discrepancy between these two results. 

In this work we identify discs with signs of grain growth 
from observations at 3, 7, 15 mm, 3-1-6 cm of 20 T Tauri 
stars in the Chamaeleon and Lupus star forming regions 
conducted using the Australia Telescope Compact Array 
ATCAq, ex tending the 3 and 7 mm work of ILommen et all 
2007I . l2010f ). An introduction to the survey and the data 
reduction process can be found in Section [5] with the re- 



The Australia Telescope Compact Array is operated by the 
Australia Telescope National Facility (ATNF) which is a division 
of CSIRO. 
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suits of the survey presented in Section [S] In Section |4] the 
dust opacity indices, the dust disc masses, and dominant 
emission mechanism at the longer wavelengths are deter- 
mined, and the potential correlation between the IR and 
mm grain growth signatures is explored, with the conclu- 
sions presented in Section [S] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample 

This survey targets 20 T Tauri stars in the Chamaeleon 
and Lupus star forming regions at a range of wavelengths 
presented in Table [1] The primary aim was to obtain 3 
and 7 mm fluxes for all sources. The so urces were se- 
lected to ove rlap with the sample observed bv lLommen et al 



ll2007l I2OIOI) and the Spitzer c2d program me (|Young et al 



I2OO5I : lAlcala et aLlliooi : iMerm et al.ll2008h . with additional 
source s wit h strong 1.3 m m fluxes from iHenning et al.l 
(|l993h and iDai et all (|2010l ) to improve the chance of de- 
tection at 3 mm. The sample includes 7 Chamaeleon I 
cloud sources, the isolated Chamaeleon source T Cha, the 
Chamaeleon II source DK Cha, 10 sources from the Lupus 1, 
3, and 4 clouds, and the isolated Lupus source RXJ1615.3- 
3255. Note that Sz 32 and WW Cha were observed simulta- 
neously at 7 mm and beyond, where Sz 32 was the primary 
target and WW Cha was in the field of view of the observa- 
tions. 

A subsets of Chamaeleon sources with the highest ex- 
pected centimetre fluxes were observed at 15 mm (6 sources) 
and 3-f 6 cm (3 sources) in order to determine the emission 
mechanisms. The expected cm fluxes were estimated by ex- 
trapolating the spectral slope from 3 to 7 mm, 03-7, using 
the 3 and 7 mm fluxes, where a = d(log F ^) / d{logiy) . No 
Lupus sources were available during the scheduled observ- 
ing time for long wavelength observations. 



2.2 Observations and data calibration 

From May 2009 to August 2011 we performed continuum 
observations with ATCA in the mm and cm bands. ATCA 
has an east-west track and northern spur with five movable 
antennas with a wavelength range from 3 mm to 30 cm, and 
one fixed antenna (CA06) at 3 km on the eastern end of 
the track with no 3 mm receivefl These observations used 
ATCA in compact hybrid configurations - observations de- 
tails and specific array configurations can be found in Ta- 
ble IBll Compact hybrid configurations were chosen to en- 
sure good u-v plane coverage and hence maximum detection 
sensitivity in relatively short integration times. 

These observations used the new Compact Array 
Broadband Backend (CABB which has a maximum band- 
width of 2 GHz (a factor 16 improvement), higher level of 
data sampling, and improved continuu m sensitivity by a fac - 
tor of four from the previous backend (|Wilson et al.ll201ll) . 



Additional information about ATCA can be found at 
http: / /www. narrabri.atnf.csiro.au/observing/ 
users_guide/html/atug.html 

^ For more information about CABB see: 
|http: / /www.narrabri.atnf.csiro.au/observing/CABB.htmI| 



For these continuum observations we used the 2048 channels 
of 1 MHz width. The observations were conducted in dual 
side band mode, with frequency pairs centred at 93+95 CHz 
(3 mm band), 43+45 GHz (7 mm band), 17+19 GHz (15 mm 
band), and simultaneous observations at 9.9-f5.5 GHz (3 
and 6 cm bands respectively). 

The gain calibrator was observed between each target 
scan of 10 minutes length at 3 and 7 mm, of 3 or 5 min- 
utes at 15 mm (see Table |D2)| and of 15 minutes at 3-1-6 cm. 
QSO B1057-79 and QSO B1600-44 were used as gain cali- 
brators for Chamaeleon and Lupus respectively, Uranus was 
used as the primary flux calibrator at 3 and 7 mm, and the 
quasar QSO B1934-638 was used as the primary flux cali- 
brator from 15 mm to 6 cm. 

The weather was generally good throughout the obser- 
vations. We reached a sensitivity of ~0.1 mjy/beam required 
to detect 90% of the sources at better than 5a for the 3 and 
7 mm band observations, and detected 4/6 sources at 15 mm 
and 1/3 at 3+6 cm. The two 7 mm non-detections had the 
lowest 1.3 and 3 mm fluxes and the highest RMS during the 
observations, and the non-detections at longer wavelengths 
had the lowest extrapolated fluxes of the targeted sample. 

The data calibration followed the standard CABB pro- 
cedure described in the ATCA user guid^l for all wave- 
lengths using th e software package MIRIAD version 1.5 
(|Sault et al.lll995l ). 



3 ANALYSIS AND RESULTS 

The flux and RMS values at 3 and 7 mm were extracted 
from the image plane using IMFIT and IMSTAT. At the longer 
wavelengths, values were extracted from the u-v plane using 
UVFIT and UVRMS due to the short 3-15 minute scans. A 
summary of the fluxes are presented in Table [2] 

The resulting source maps were created by Fourier- 
transforming the visibility data using a speciflc weighting 
(here "natural"), then the CLEAN and RESTORE algorithms. 
For sources observed over several epochs, data from each 
day was calibrated separately and combined with INVERT. 



3.1 Which sources are resolved at 3 and 7 mm? 

We determined if a source was extended by looking at the 
visibilities. As the resolution of the array increases, the flux 
will decrease for a resolved source and remain constant for 
an unresolved source. Plots of the u-v distance as a function 
of amplitude created using the UVAMP task in MIRIAD are 
presented in Fig. [2] for 3 mm and Fig. |3] for 7 mm. These 
plots suggest DK Cha at 3 and 7 mm and CS Cha at 3 mm 
have extended emission, while marginally extended emission 
is observed at 7 mm for CR Cha, CS Cha, SZ Cha and 
RXJ1615. 3-3255. For these marginally resolved sources, we 
will take the conservative approach and assume the emission 
were not resolved. 

RXJ1615. 3-3255 was resolved at 1.3 mm with the SMA, 
with an estimated Gaussian siz e of 1.53 arcsec (~28 AU as- 
suming a distance of 184 pc) (|Lommen et al.|[201(J l. Since 



* http://www.narrabri.atnf.csiro.au/obscrving/uscrs_guide/ 
html/atug.html 
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Table 1. List of sources observed with ATCA in this survey. 



Source 


RA 


DEC 


Spectral Type 


Cloud 


Distances 


Toff 


Wavelengths 


Comments 




(J2000) 


{J2000) 






(pc) 


(K) 


(mm) 





Chamaelcon 



1 


SY Cha 


10 56 30.5 


-77 11 39.3 


M0.5 


Cha I 


160 


3778 


7 




2 


CR Cha 


10 59 06.9 


-77 01 39.7 


K0-K2 


Cha I 


160 


4900 


7,15 




3 


CS Cha 


11 02 24.9 


-77 33 35.9 


K4 


Cha I 


160 


4205 


3,7,15 


Binary" 


4 


DI Cha 


11 07 21.6 


-77 38 12.0 


G1-G2 


Cha I 


160 


5860 


3,7,15 


Binary" 


5 


T Cha 


11 57 13.6 


-79 21 31.7 


G2 


Isolated 


100 


5600 


7,15,30,60 




6 


Glass I 


11 08 15.1 


-77 33 59.0 


K4 


Cha I 


160 


5630 


3,7 


Binary" 


7 


SZ Cha 


10 58 16.7 


-77 17 17.1 


KO G 


Cha I 


160 


5250 


7 




8 


Sz 32 


11 09 53.4 


-76 34 25.5 


K4.7 


Cha I 


160 


4350 


3,7,15,30,60 


HH jet" 


9 


DK Cha 


12 53 17.2 


-77 07 10.7 


FO D 


Cha II 


178 


7200 


3,7,15,30,60 


Outflow" 



Lupus 



10 


IK Lup 


15 39 27.8 


-34 46 17.2 


K7 D 


Lupus 1 


150 


3802 


7 


Binary'' 


11 


Sz 66 


15 39 28.3 


-34 46 18.0 


M2 D 


Lupus 1 


150 


3350 


7 


Binary*" 


12 


HT Lup 


15 45 12.9 


-34 17 30.8 


K2 


Lupus 1 


150 


4898 


7 


Triple system, HH jet'' 


13 


GW Lup 


15 46 44.7 


-34 30 36.0 


M2 


Lupus 1 


150 


3499 


7 




14 


GQ Lup 


15 49 12.1 


-35 39 03.9 


K7V D 


Lupus 1 


150 


3899 


3,7 


Planet'' 


15 


RY Lup 


15 59 28.4 


-40 21 51.2 


GOV: C 


Lupus 3 


200 


4592 


7 




16 


HK Lup 


16 08 22.5 


-39 04 46.3 


K5 


Lupus 3 


200 


4350 


7 




17 


Sz 111 


16 08 54.7 


-39 37 43.1 


Ml. 5 


Lupus 3 


200 


3573 


7 


Cold disc 


18 


EX Lup 


16 03 05.5 


-40 18 25.3 


MO 


Lupus 3 


200 


3802 


3,7 




19 


MY Lup 


16 00 44.6 


-41 55 29.6 




Lupus 4 


165 


5248 


7 




20 


RXJ1615. 3-3255 


16 15 20.2 


-32 55 05.0 


K5 


Isolated 


184 


4590 


7 





" See Section [AT] for further details on Chamaeleon sources. 
'' See Section IA2I for furt her details on Lupus so urces. 
References: Chamaeleon: jHughes et"al] lll99ll ') : IWhittet et al.l lll997l} : 



2007); 



19971); 



.Spezzi et al.| (120081) . Lupus: IComeronI ||2008 
Feldt et al.l l|l998h 



van den Ancker et al.l lll998l): iLuhmanI 


2004^: iBrown et al. 


Mortier et alj ll201lh: Ivan Kempen et al. 


ll2010l);IChen et al. 



the H168 array has a resolution of 1753 x 763 AU (beam 
size of 9.5x4.1 arcsec) and the map in Fig. |4a]does not ap- 
pear to have extended emission, it is unlikely the emission 
was resolved here at 7 mm. 

The visibilities for Sz 32 and MY Lup presented in 
Fig.[3]suggest the presence of a contribu ting source near each 
target ed source location. Following the lFomalont fc Wrightl 
1 19741 ) approach of using the visibility function to determine 
the distance between double sources, we find that the Sz 32 
and MY Lup visibilities are consistent with the location of 
WW Cha at ~ 35" to the south-west of Sz 32 and an uniden- 
tified source at ~ 50" to the north of MY Lup - maps are 
presented in Fig.U] The smaller field of view at 3 mm (30". 2 
versus 65". 4 at 7 mm) did not allow for either WW Cha or 
the unident ified source to be dete cted at 3 mm. This is con- 
sistent with iLommen et al.l (|2009l ) , who previously observed 
and resolved WW Cha at 3 and 7 mm and only detected Sz 
32 in the field of their 7 mm observations. 

Although looking at the visibility data is commonly 
used in literature to dete rmine source extension (e.g. 
[Andrews fc Williamsl l2007bl ) . this approach could suffer 
from phase decorrelation. Phase decorrelation degrades the 
signal as the baseline increases, causing a similar decrease in 
amplitude on the longer baselines. Given the generally good 
weather conditions, the use of compact array configurations. 



and reasonable RMS values from the ATCA seeing monitoiQ, 
phase decorrelation was not substantial in our observations. 

To ensure that phase decorrrelation was not effecting 
our visibilities amplitudes, we used three other methods to 
investigate source extension at 3 and 7 mm. If a source 
is extended, a point source fit will exclude some emission, 
whilst a Gaussian fit should include most of the emission. 
A source could be considered extended when (1) the point 
source fit is less than the Gaussian fit (after accounting for 
uncertainties), (2) the Gaussian fit is greater than the point 
source fit plus n times the RMS, w here n is a predetermined 
value (|Lommen et al.1 (|2007l . |2009| ) used n = 2 for pre-CABB 
ATCA data, and assuming a factor of four improvement 
with CABB, here we use n = 8) and (3) when the syn- 
thesised beam is smaller than the Gaussian size obtained 
from the Gaussian fit. These three alternative approaches 
were considered (see Table IC1|) and the results are consis- 
tent with the sources considered resolved when using the 
visibility data (CS Cha and DK Cha), while marginally re- 
solved sources were considered unresolved with these three 
alternative approaches. 

Thus from this analysis we have determined that the 
emission of CS Cha was resolved at 3 mm, and the emission 
of DK Cha was resolved at 3 and 7 mm. In Table [5] sources 

^ The ATCA seeing monitor is a fixed interferometer tracking the 
30.48 GHz beac on on a geostationary co mmunication satellite at 
elevation of 60° jMiddelbere et al.llioO^'l . 
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in boldface were r esolved at 3 mm with ATCA in this work 
or in the literature. iLommen et all (|2007l . [2"oiCf ) resolved SZ 
Cha, HT Lup, RY Lup and Sz 111 at 3 mm with ATCA. 



3.2 Source fluxes 

A summary of all the continuum fluxes obtained in this sur- 
vey are presented in Table (2] The flux values were obtained 
by combining both sidebands for all available days. The Ta- 
ble is suppl emented with 1 . 2 mm fluxes for Chamaeleon 
sources from iHenning etHI (liggsl) wi th S EST, and for Lu - 



pus sources from Lommen et al.l ( 2010[) and lDai et al.l (|201Cll ) 



with the S ubmillimetre A r ray (SMAn , as well as some 3 mm 
data from lLommen et al.l (|2007l . 12010 ) with ATCA - see Ta- 
ble [T] for a list of sources observed and wavelengths used in 
this survey. Four of our Chamaeleon s ources have 870 ^m 
LABOCaQ data jBelloche et al.ll2011bl '). however for consis- 
tency we will use the SEST flux data for our analysis of the 
Chamaeleon sources. 

All seven sources observed at 3 mm in this survey were 
detected, and are consistent with the pre vious upper lim- 
its reported in ILommen et al.l (|2007l . l2010l l . We present for 
the first time observations of the sources at 7 mm (for excep- 
tions see Table IdTI) and beyond, 18/20 sources at 7 mm, 4/6 
sources at 15 mm and 1/3 sources at 3-1-6 cm were detected. 
A 3(T upper limit is provided for all non-detections. 

If the source was considered resolved in Section 13.11 or 
in the literature, the Gaussian fit for the flux is presented 
in Table [2] along with the flux uncertainties in the fits, the 
RMS values, and the beam size for observations made in 
this survey. The primary flux calibration uncertainties are 
not included in Table[2] This uncertainty is '-20% at 1.2 mm 
for SEST data and -10% for SMA data. For the ATCA data 
this uncertainty is ~30% at 3 mm and ~10% at 7-15 mm and 
3-1-6 cm bands. 

The beam size of 4". 85 x 3". 53 at 7 mm was sufficient 
to resolve the binary pair IK Lup and Sz 66 with a sepa- 
ration of 6". 4 (|Reipurth fc Zinneckeij|l993l ). However, Sz 66 
was not detected - see Fig. I4bl Thus the 7 mm fluxes in 
Tables ICll are most l ikely only frorn IK L up. Similar re- 
sults were reported bv I Lommen et al] (|201(j ) for the 1.2 and 
3.2 mm observations, with beam sizes of 5".0 x 2".l and 2" 
respectively. In both cases the beam size was sufficient to 
resolve the binary, su ggesting the 2.2 mjy flux reported by 
ILommen et"all (j20ld ) at 3.2 mm for Sz 66 is likely an over- 
estimate. If Sz 66 was not detected at 3.2 mm, a 3a" upper 
limit of 1.2 mJy would be more appropriate. Sz 66 was de- 
tected with Sp i tzer a nd has a 10 /im silicate feature - see 
ILommen et al.l (|2010l l. However, the lack of cold dust from 
1.2 to 7 mm would suggest there is no circumbinary disc, 
and if Sz 66 has a cold dust disc its mass is too low to have 
been detected by our observations. 



^ The Submillimetre Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and funded by the Smithso- 
nian Institution and the Academia Sinica. 

LABOCA: LArgc BOIometer CAmera - service at the 12-m 
APEX telescope 



6 Ubach et al. 



20 40 60 80 100 120 
U-V distance (klom) 



2,5 

2.0 
1 .5 
1 ,0 
0,b 
0.0 



20 40 60 SO 100 -20 
U-V distonce (klarn) 



2.5 ~ 
2.0 : 

1 .5 : 
1 .0 : 
0.5 : 

0.0 _ 



20 40 60 SO 100 120 
U-V distance (klam) 




20 40 60 80 100 120 
U-V distance (klam) 




20 40 60 80 100 '20 
U-V distance (klam) 



OQ Lup 



20 40 60 

U-V distance (klam) 



2.5 
„ 2.0 
I 1,5 

■S 1 ,0 

< 

0,5 
0,0 



EX Lup 



40 60 SO 100 -20 
U-V distdnce (klam) 



Figure 2. Visibility amplitude versus baseline length (or u-v distance) for Chamaeleon and Lupus sources at 3 mm, with the Icr statistical 
error bars for each bin. 
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Figure 3. Visibility amplitude versus baseline length (or u-v distance) for Cliamaeleon sources at 7 mm, with the Icr statistical error 
bars for each bin. 



8 Ubach et al. 



10 20 30 40 50 
U-V distance (klom) 



HT L.ip 



10 20 30 
U-V distonce (klarn) 



1 .0 
O.B 
0.6 
0.4 
0.2 
0.0 



OQ Lup 



10 20 30 40 
U-V distonce (kldm) 




1 20 30 

U-V distance (klam) 




10 20 30 40 bO 

U-V distonce (klan) 



10 20 30 

U-V distance (klam) 



10 20 30 40 

U-V distonce (klam) 



■s 2 ; 

Q- 

E 
< 

1 f 





ivir Lup 



10 20 30 40 
U-V distdnce (klam) 



1 .2 
1 .0 
0.8 
0.6 
0.4 
0.2 



10 20 30 40 50 

U-V distonce (klam) 



EX Lup 



20 40 50 

U-V distance (klam) 



Figure 3. Visibility amplitude versus baseline length (or u-v distance) for Lupus sources at 7 mm, with the Icr statistical error bars 
each bin. 
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Figure 4. The cleaned maps using natural weighting for RXJ1615. 3-3255 at 7 mm, IK Lup at 7 mm, ATCA160043-415442 at 3 mm and 
MY Lup at 7 mm. (A) RXJ1615. 3-3255 at 7 mm with contours at -3,3,5 times the image RMS of 0.2 mjy/beam. (B) IK Lup (Sz 65) 
at 7 mm with contours at -3,3,6,9,12 times the image RMS of 0.8 mJy/beam. Although the synthesised beam size of 4". 9 x 3". 5 was 
sufficient to resolved the binary with a separation of 6". 4, Sz 66 was not detected. (C) ATCA160043-415442 at 3 mm (with a cross in 
the position of MY Lup) with contours at -3,3,12,21,31 times the image RMS of 0.3 mJy/beam. (D) Beam flux corrected clean map of 
MY Lup at 7 mm with ATCA160043-415442 in the north with contours at -3,3,5,11,22,45 times the RMS of 0.3 mjy/bcam. 



Table 2. Survey results. (1) Source name. (2) 1.2 mm continuum flux with image RMS in parentliesis. (3) 3 mm flux with image RMS in parenthesis. (4) Beam size at 3 mm. (5) 7 mm 
flux with image RMS in parenthesis. (6) Beam size a t 7 mm . (7), (8), (9) F luxes at 15 mm, 3 and 6 cm with u-v plane RMS in parenthesis. flO) References : For 1 .2 mm continuum: 
1-Hen ning et al. (1993) with SEST, 2 4Lommen et alj ll2010h - with SMA, 3- |Dai et all | |2010|'| with SMA. For 3 mm continuum with ATCA: 4 4Lommen et al.l ||2007| '), 5- |Lommen et alj 
1I2OIOI '). 6-this work. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Source 


Fl.2mm (RMS) 


Fsmm (RMS) 


Beam Size 


F7mm (RMS) 


Beam Size 


FlSmm (RMS) 


Fscm (RMS) 


Fscm (RMS) 


References 




mjy (mjy/beam) 


mJy (mjy/bcam) 


arcsec 


mJy (mJy/beam) 


arcsec 


mJy (mJy/beam) 


mJy (mJy/beam) 


mJy (mJy/beam) 




Chamaclcon 


SY Cha 


<172.0 


<4.8 




<0.9 (0.3) 










1,4,6 


CR Cha 


124.9 (24.2) 


6.2 (1.5) 


2.5x2.1 


1.7 (0.1) 


5.0x4.3 


0.3 (0.1) 






1,4,6,6 


CS Cha 


128.4 (45.6) 


9.4 (0.6)°* 


13.2x1.4 


1.5 (0.1) 


5.0x4.3 


<0.3 (0.1) 






1,6,6,6 


DI Cha 


38.0 (11.4) 


2.3 (0.3)*" 


4.6x1.8 


0.9 (0.1) 


4.8x4.2 


<0.3 (0.1) 






1,6,6,6 


T Cha 


105.0 (17.7) 


6.4 (1.0) 


2.5x2.4 


3.0 (0.1) 


5.1x4.1 


0.3 (0.1) 


<0.3 (0.1) 


0.3 (0.1) 


1,4,6,6,6,6 


Glass I 


69.9 (22.4) 


4.4 (0.1)* 


2.4x1.8 


0.7 (0.1) 


4.9x4.1 








1,6,6 


SZ Cha 


77.5 (20.3) 


5.8 (0.5)° 


2.3x2.1 


0.7 (0.1) 


5.0x4.1 








1,5,6 


Sz 32 


93.1 (20.8) 


3.1 (0.2)* 


2.3x1.8 


1.2 (0.1) 


5.0x4.0 


0.4 (0.1) 


<0.3 (0.1) 


<0.3 (0.1) 


1,6,6,6,6,6 


DK Cha 


680.0 (22.0) 


49.8 (1.3)°* 


6.4x1.3 


6.6 (0.6) 


47.3x2.5 


1.6 (0.1) 


0.8 (0.1) 


<0.3 (0.1) 


1,6,6,6,6,6 


Lupus 


IK Lup 


28.0 (2.8) 


3.4 (0.4) 


2.1x1.7 


0.9 (0.1) 


4.5x3.5 








2,5,6 


Sz 66 


<8.0 (2.8) 


2.2 (0.4)'' 


2.1x1.7 


<0.3 (0.1) 


4.5x3.5 








2,5,6 


HT Lup 


73.0 (4.0) 


12.0 (1.1)° 


2.4x1.7 


3.4 (0.1) 


4.4x3.5 








4,4,6 


GQ Lup 


25.0 (3.0) 


3.6 (0.3) 


2.8x2.0 


0.6 (0.1) 


6.7x4.4 








3,6,6 


GW Lup 


64.0 (3.7) 


8.5 (1.9) 


5.3x1.7 


0.6 (0.1) 


4.8x3.5 








4,4,6 


RY Lup 


89.0 (4.9)° 


2.8 (0.7)° 


2.2x1.9 


1.0 (0.1) 


4.7x3.5 








2,5,6 


HK Lup 


101.0 (3.9)° 


7.3 (2.1) 




1.0 (0.1) 


5.2x3.5 








4,4,6 


Sz 111 


52.5 (4.2)° 


5.7 (0.7)° 


2.1x1.9 


0.5 (0.1) 


5.2x3.5 








5,5,6 


EX Lup 


21.3 (4.0)° 


2.0 (0.3) 


2.9x1.9 


1.4 (0.1) 


4.9x3.6 








2,6,6 


MY Lup 


66.1 (3.4)° 


8.7 (0.4) 


2.0x1.7 


1.1 (0.1) 


4.8x3.4 








2,5,6 


RXJ1615. 3-3255 


169.1 (3.9)° 


6.7 (0.6) 


2.0x1.7 


0.8 (0.2) 


9.5x4.1 








2,5,6 



CI 



A 3cr upper limit is provided for non-detections. 
Sources in boldface were resolved at 3 mm with ATCA. 
" Flux at 93 GHz, as 95 GHz was not detected. 
' Potentially an over estimate - see Section 13.21 

* 3 mm continuum fluxes from this work. All 7 mm fluxes are from this work. 

° Gaussian flux flt for sources resolved at 1 mm with SMA and 3 mm with ATCA. 
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3.3 Millimetre spectral slopes 

The mm spectral slopes a and 03-7 presented in Table [3] 
were determined using the 1, 3 and 7 mm band fluxes from 
Table [2l The uncertainties in a and 03-7 were calculated 
through propagation of the flux uncertainties presented in 
Table ICll Histogram of the distribution of a values is pre- 
sented in Fig. [5l The values of a range between 2 to 4. 
The grain size distribution in disc models evolves and can 
depart substantially from the n(a) oc a~^'^ size distribu- 
tion of the interstell ar medium (e.g, IWeidenschillind 1 19971 : 
iBirnstiel et~all I2OIII ). Since fully modelling the grain size 
distribution for all 20 sources in our survey is beyond the 
scope of this paper, for this analysis we will assume the dust 
emissivity can be represented by a single power law from mm 
through to cm wavelengths. Thus if only thermal dust emis- 
sion is contributing to the flux from 1 to 7 mm the spectral 
slope would remain constant, whilst an inequality between 
a and Q!3-7 would imply a break in the spectral slope, sug- 
gestive of the presence of other emission mechanisms such 
as thermal free-free and non-thermal emissions contributing 
to the flux at the longer wavelengths. 

For this analysis we will consider a and 03-7 to be con- 
sistent (no break at 7 mm) when Q!3-7 ± Aa^-r = a± Aa and 
inconsistent (a break at 7 mm) when 03-7 ± Aq:3_7 < a±Aa. 
The uncertainty Aa is 0.4 for all sources and Aq3_7 is given 
in Table El 

The histogram a values show that majority of Lupus 
sources have a value of q < 3, while majority of Chamaeleon 
sources have a value a > 3 - see Fig. (5] From Table [3] we 
determined that 11 sources do not have a break in their 
spectral slope at 7 mm, while 7 sources do have a break and 
two are undetermined. 

We plot the flux as a function of wavelength in Figs. [6] 
and [7] for Chamaeleon and Lupus respectively, to better 
study the spectral slope. The solid line represents the spec- 
tral slope a presented in Table [3] and the dashed-dot lines 
are the upper and lower limits for a using the flux fit uncer- 
tainties and assuming the standard primary flux calibration 
uncertainties (20% for 1.2 mm, 30% for 3 mm, and 10% for 
7, 15 mm and 3, 6 cm). Thus the shaded region represents 
the flux range expected when thermal dust emission domi- 
nants. T he dashed line represe nting free-free emission with 
Qfi = 0.6 (|Panagia fc Fellll975l ) anchored at the 15 mm data 
point is included for sources with cm data. Full SEDs can 
be found in Figs. [CTI and [C2l 

CR Cha, DI Cha, T Cha, Sz 32, RY Lup, EX Lup and 
RXJ1615. 3-3255 all have 7 mm fluxes in excess of purely 
thermal dust emission consistent with Table [3]- see Figs. [S] 
and[7l Note that if we were to use t he 870 fim fluxes obtained 
with LABOCA bv lBelloche et all (|2011al ') for the three avail- 
able Chamaeleon sources instead of the SEST data, the a 
values would change to 2.3±0.2 for CS Cha, 2.9±0.2 for Glass 
I and 3.1±0.4 for SZ Cha, which would not change the in- 
terpretation of our analysis of these sources. 

For those sources with observations at 15 mm and 
longer wavelengths, we can see that all six Chamaeleon 
sources observed at 15 mm have excess emission above that 
expected from thermal dust. At 3-1-6 cm, DK Cha, T Cha 
and Sz 32 the fluxes and upper limits suggest the presence 
of free- free emission from an ionised wind and/or chromo- 



Table 3. Spectral slopes from 1-3 mm and 3-7 mm. (1) Source 
name. (2) Spectral slope from 1-3 mm, a (with an uncertainty of 
0.4). (3) Spectral slope from 3-7 mm, Q3-7. (4) Uncertainties in 
Cfi-7- (5) Indication of whether there is a break in the spectral 
slopes at 7 mm. 



Sources 


a 


03-7 


A«3-7 


Break 


Chamaeleon 


SY Cha* 


>3.9 


>2.1 


0.5 




CR Cha 


3.4 


1.7 


0.4 


Y 


CS Cha 


2.9 


2.5 


0.3 


N 


DI Cha 


3.2 


1.2 


0.3 


Y 


T Cha 


3.1 


1.1 


0.3 


Y 


Glass I 


3.1 


2.5 


0.4 


N 


SZ Cha 


2.9 


2.9 


0.5 


N 


Sz 32 


3.8 


1.3 


0.3 


Y 


DK Cha 


2.9 


2.7 


0.2 


N 


Lupus 


IK Lup 


2.4 


1.7 


0.4 


N 


Sz 66* 


>1.4 


>3.3 


0.3 




HT Lup 


2.0 


1.7 


0.4 


N 


GQ Lup 


2.2 


2.4 


0.5 


N 


GW Lup 


2.3 


3.6 


0.4 


N 


RY Lup 


3.9 


1.4 


0.4 


Y 


HK Lup 


2.9 


2.7 


0.4 


N 


Sz 111 


2.5 


3.3 


0.5 


N 


EX Lup 


2.6 


0.5 


0.2 


Y 


MY Lup 


2.3 


2.8 


0.5 


N 


RXJ1615. 3-3255 


3.6 


2.8 


0.5 


Y 



* SY Cha and Sz 66 only have upper limits at 1.2 and 3 mm. 
The sources in bold were resolved at 3 mm. 



8 







■ 

□ 


Cham 
jupus 






> 2.5 : 


3.5 4 



a 



Figure 5. Histogram of the spectral slopes a. SY Cha and Sz 66 
spectral slopes are unconstrained due non-detection at 1.2 and 
3 mm, thus both are excluded from the analyses. 

spheric emission, which will be discussed further in the next 
Section. 

3.4 Emission mechanisms at longer wavelengths 

Temporal flux monitoring helps to determine the presence 
of contributing emission mechanisms. For our sample, a sub- 
set of 4 Lupus and 2 Chamaeleon sources which have 7 mm 
fluxes obtained approximately one year prior to our survey 
- see Table [DTI allowing us to determine if free- free emission 
is contributing to the 7 mm fluxes for these sources. We also 
observed a sub-sample of 6 Chamaeleon sources at 15 mm 
and 3 of these sources at 3-1-6 cm for one epoch (see Table[2]). 
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Figure 6. Millimetre flux versus wavelength for Chamaeleon sources. The solid line represents the spectral slope a and the dashed-dot 
lines are the upper and lower limits for a using the flux fit uncertainties and assuming the standard primary flux calibration uncertainties 
(20% for 1.2 mm, 30% for 3 mm, and 10% for 7 and 15 mm and 3+6 cm). The dashed line represents the free-free emission with Off = 0.6 
anchored at the 15 mm data point, is included for sources with cm data. Fluxes values from Table [2] Note that the error bars include 
the flux uncertainty from the fitting routine and the primary flux calibration uncertainties. 



We can investigate flux variations in these sources by deter- 
mining a point source fit in the u-v plane for each scan length 
at these wavelengths allowing us to determine the contribut- 
ing emission mechanisms 15 mm wavelength and beyond. 
The results from the long wavelength temporal monitoring 
are presented in Figs. [8] and [9] with the maximum and min- 
imum recorded flux values at 15 mm and 3+6 cm presented 
in Tables ID2I and ID3I respectively. 

The liter a ture 7 mm data were obtained from 
iLommen et al. ] (|2009l . l2010l ). Their observations took place 
in 2008 with the ATCA using the pre-CABB system. Their 
observing frequency pairs at different epochs do not always 
match our central frequencies. However small changes in the 
observing frequency (42.8 to 43.0 GHz and 45.1 to 45.0 GHz) 
between different epochs would not result in a change in flux 
by more than the statistical uncertainty and thus should not 
affect our conclusions. 

For our subsample at 7 mm, we can also investigate the 
potential day-to-day flux variability by obtaining a point 
source fit fiux for the central frequency of the dual side- 



band for each day the source was observed, resulting in in- 
tegration times between 10 to 70 minutes depending on the 
source (see Table [Blj) . In general, sources observed on con- 
secutive days do not have the same integration times, and 
since Ghamaeleon never sets at ATCA, more time was spent 
on the Ghamaeleon than the Lupus sources. 

For this analysis we define "short temporal monitoring" 
as scans separated by less than a day, and "long temporal 
monitoring" as scans separated by a day or more. Thermal 
dust emission is considered dominant when no intra-epoch 
variability is observed. Thermal free-free emission is likely 
present when the fiux varies during long temporal monitor- 
ing by a factor of 20-40%. Non-thermal emission will be 
considered present when the flux varies by a factor of 2 or 
more on timescales of minutes to hours. For this analysis 
we take into account the flux flt uncertainty and the pri- 
mary flux calibrator uncertainty for GABB data. We only 
claim variability if the same behaviour is observed in both 
frequency sidebands. 
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Figure 7. Same as Fig. [G] but for Lupus sources. 

3.4-2 Chamaeleon sources 



Four Lupus sources (RY Lup, MY Lup, RXJ1615. 3-3255 and 
Sz 111) were observed three times at 7 mm, once in 2008 and 
twice in 2009 on consecutive days - fluxes are presented in 
Table iDl] RY Lup, MY Lup and Sz 111 show no evidence 
of flux variability over daily or yearly timescales, indicating 
that thermal dust emission is dominant, while RXJ1615.3- 
3255 is consistent with a presence of thermal free-free emis- 
sion. The a values for all sources are consisten t with the 
upper limits calculated bv lLommen et all (|2010l ). 



Temporal flux monitoring at 7 mm is available for CS Cha 
and Sz 32 - see Table IDll The 15 mm temporal monitor- 
ing for 6 Chamaeleon sources are presented in Fig. [5] with 
the maximum and minimum recorded flux presented in Ta- 
ble ID2I The 3-1-6 cm temporal monitoring for 3 of those 6 
Chamaeleon sources are shown in Fig. [5] with the recorded 
fluxes in Table [D3l 



The temporal monitoring of CR Cha and DI Cha sug- 
gest that the excess emission seen in Fig. [6] at 15 mm is not 
from a fast varying chromospheric emission but most likely 
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from thermal free-free emission, consistent with the spectral 
slope of 0.6 seen in Fig. |6l 

CS Cha was observed at 7 mm three times in 2008 with 
extended array configurations (1.5 and 6 km), and twice in 
2009 on consecutive days with a compact array configura- 
tion. There is no evidence of flux variability at 6.7 mm on 
consecutive days in 2009 (see Table IDll) . while the three 
fluxes obtained in 2008 are inconsistent with each other and 
with the fluxes obtained in 2009. The 2008 and 2009 changes 
are consistent with the fact that CS Cha was resolved with 
the extended arrays in 2008 and unresolved with the com- 
pact array in 2009, which likely caused an underestimate 
of the flux when CS Cha was resolved. The variability in 
the 7 mm flux observed in 2009 is consistent with thermal 
free-free emission, and the likely cause of the excess emission 
observed in Fig. [S] 

The temporal monitoring of DK Cha at 15 mm is consis- 
tent with the presence of some excess emission from thermal 
free- free emission - see Fig. [8] and Table [D2l - and contami- 
nation from chromospheric emission was detected at 3-f 6 cm 

- see Fig. and Table lD3l These results are consistent with 
the excess emission at 15 mm and 3-1-6 cm seen in Fig. [6] 

Sz 32 was observed once in 2008 and twice in 2009 on 
consecutive days at 7 mm. The level of flux variability on 
consecutive days and between the 2008 and 2009 (see Ta- 
ble ID1[I indicate thermal free-free emission is contributing 
to the 7 mm flux. 

T Cha and Sz 32 have excess emission at 15 mm likely 
due to contamination from thermal free-free emission from 
an ionised wind, and excess emission at 6 cm from chromo- 
spheric emission — see Tables ID2I ID3I and Figs. [8] and (9] 
These results are consistent with the excess emission seen in 
Fig. El This is also consistent with Sz 32 being the suggested 
source that drives the short E ast- West jet, HH 914 se en at 
a distance of 0.4" from Sz 32 (|Wang fc Henningj|2006l ). 

3.5 Additional Sources 

WW Cha was detected in the field of view of the Sz 32 
observations at 7, 15 mm and 3-1-6 cm. The continu um fluxes 
obtained are consistent with lLommen et al.l (|2009l ). 

A second source was detected at 7 mm to the north- 
west of MY Lup - see Fig. I4dl for beam flux corrected 
clean map. This is an unknown source which we identify 
by ATCA160043-415442 (instrument + J2000 location), first 
observed in May 2009 at RA 16:00:43, Dec -41:54:42 with an 
estimated Gaussian flux of 25.0 ± 5.5 mjy after a beam flux 
correction. Since the source was at the edge of the beam, 
observations were conducted at 3 mm to confirm the detec- 
tion and investigate the spectral slope. In August 2010, the 
source was detected in the same location at 3 mm with a 
point source fit of 9.0 ± 0.4 mJy and RMS of 0.3 Jy/beam 

- see Fig. [Ic] There are no known 1.2 mm or cm band de- 
tections in the literature. The mm spectral slope from 3 to 
7 mm, 03-7 = -1.4 ± 0.8, suggests non-thermal emission and 
the source is likely a background radio galaxy. 

A range of extra-galactic sources were detected during 
the 3+6 cm observations, including: SUMS J115749- 792819, 
SUMS J11558-792403, SUMS J1201169-790726, G 2MASX 
J11574287-7923537, PMN JllO-7618, PMN J1113-7625, 
SUMSS J110856-763624, CXO J110841.4-763451, SUMS 
J125216-770903 and SUMS J125225-770345. 



4 DISCUSSION 

In this section we aim to identify the protoplanetary discs 
with signs of mm and cm-sized grains. Thus far we have 
determined the sources with extended emission (i.e. which 
we assume to be from resolved discs) in the 3 mm band, 
breaking the degeneracy between radially extended optically 
thin discs and compact optically thick discs with shallow 
spectral slopes. 

Analysis of the continuum fiuxes in the mm and cm 
bands and the resulting spectral slopes allow us to identify 
sources with a break in spectral slope at 7 mm, indicating 
contributions to the emission at 7 mm and beyond from 
other emission mechanisms other than thermal dust. Here 
we calculate the dust opacity index at 3 mm in order to 
estimate the maximum grain sizes at this wavelength, and 
then estimate the dust disc masses at 3 mm. Finally, we 
explore a previously claimed tentative correlation between 
grain growth signatures in the IR and mm regimes using this 
data set. 

4.1 Dust opacity index at 3 mm 

The 1 and 3 mm band emission are in the Rayleigh- Jeans 
limit and is expected to be mostly optically thin. The spec- 
tral slope a can be used to estimate the dust opacity index 
p. Assuming the emission is optically thin, the dust opacity 
index can be written as P ~ a-2 and an opacity index /3 < 1 
indicates grain growth up to mm sizes (|Drainell200a ). Ta- 
ble 3] presents the /3 values determined assuming optically 
thin emission for all the sources in our sample. We found 
that 2 of the 6 resolved sources and 5 of the 12 unresolved 
sources have /3 < 1, suggestive of grains up to mm sizes - see 
columns two and three of Table |4l Due to the uncertainties 
in the /3 values (determined through the propagation of the 
a uncertainties) the symbol "?" is used when 0.7 < /3 < 1.4 
in Table g] 

This approximation of jS does not account for the con- 
tributi on from optically thick emission. Taking this into ac- 
count, iBeckwith erall (|l990l ') showed that 

P . (a-2)(l + A), (1) 

where A is the ratio of optically thin to thick emission given 
by 

A « -px[(2-q) ln(l-p/2)r]-\ (2) 

where q and p are the power-law exponents of the disc tem- 
perature and surface density profiles respectively, and f is 
the average disc opacity at the specific frequency. The log- 
ari thmic dependence on A is only valid when 2-p-q = 0- 
see IBeckwith erall (|l990l 'l for details. 

We determined the temperature profile exponent q for 
our sources using IRAS 60 and 100 fim fiuxes obtained from 
the Infrared Science Archival using q = 2/(3 - aiajfl, and 
found 0.2 < q < 0.7 with an uncertainty of 0.4, with the 
exception of HT Lup and GQ Lup with g > 2 — see Table ID 

* Housed at the Infrared Science Archive (IRSA) at the Infrared 
Processing and Analysis Centre (IPAC), California Institute of 
Technology. 

^ Whore a m = d(log Fty)/d(log u) is the IR spectral slope for 
A < 100 fim llBeckwith et aI.lll99ol : iRodmann et al.ll2006l) 
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Figure 8. Temporal monitoring of 15 mm band continuum flux at 17 and 19 GHz of 6 Chamaclcon sources. The scan lengths were 
between 3 and 5 minutes with total integration times given in Table IBll The fluxes are represented by triangles with corresponding 
uncertainties, and the red squares represent 3cr detection values. Note that when a source is not detected (fluxes below the red squares), 
UVFIT can sometimes estimates a negative flux. These are included for completeness, however the value itself is incorrect. 



For a fi at spectrum disc a val ue of g = 0.5 and p = 1.5 
is expected (|BecltwitIi et al.lll99Cf ). an d jiigii angular res- 
olution observat ions suggest p < 1.5 (|Wilner et al.l l2000l : 
iTesti et al-lbOOd). If w e adopt p = 1.5, a f = 0.02 at 3.3 mm 
I Lommen et al.l 2007h and use the q values in Column 4 of 



Table S] with uncertainties such that the q values are con- 
sistent with 0.5, A ~ 0.2 (with exception of HT Lup and GQ 
Lup where A approximation in Eq. [2]is invalid). 

There are three sources (HT Lup, HK Lup, and MY 
Lup) for which Spitzer 24, 60, 70, and 100 ^m fluxes ex- 
ist. Using the MlP d'^ fluxes at 24 and 70 fim obtained by 
[m erin et al. ] (|2008D , we re-calculated q and found values of 
0.4, 0.9, and 0.6 for HT Lup, HK Lup, and MY Lup respec- 
tively. The significant decrease of q from 2.9 to 0.4 for HT 
Lup is well above the 0.4 uncertainty level, while the HK 
Lup and MY Lup increases are within the uncertainty. The 
increase in q for HK Lup and MY Lup leads to no signif- 
icant change in the A and /3 values, while the decrease in 
HT Lup to Q = 0.4 allows for the use of the A approximation 
(obtaining a A = 0.2). 

Since knowledge of the inner and outer radius of the 
disc is needed to determine an exact value for A (see eq. 20 



Multiband Imaging Photometer for Spitzer is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space Adminis- 
tration. 



of iBeckwith et al.|[l990h and taking A = 0.2 changes /3 by 
0.3 or less which does not change the interpretation of our 
analysis (with exception of CR Cha which increases to 1.6 
and thus would be considered to have ISM sized grains); for 
further analysis we will use the values for /3 calculated in 
Tabled 



4.2 Dust disc masses 

The dust disc mass can be determined via 
F 

Md.sr = (3) 

using the distances D presented in Table [T] the 3 mm fluxes 
Fv from Table [21 the brightness, i3i;(rdust), for a dust tem- 
perature Tdust (assumed to be 25 K) given by the Planck 
function, and the dust opacity, k„ at tj = 94 GHz, evaluated 
using the P values in Table [4] via the IBeckwith et"al] (|199G| ) 
equation 



10(i;/10^^ Hz) 



J 2 -1 

cm g . 



(4) 



The resulting Ksmm and Mduat values are presented in 
Table |4l Note that Ksmm varies from the canonical 
0.9 cm^g''^ used in the past works (e.g., Beckwith et al.l 
1 19901 : 1 Andrews fc WiUiams|[2007al : iRicci et al.ll2010bD . 

The dust disc masses which range from lO^'^-lO"'' M© 
are similar to those found in other star forming regions (e.g.. 
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Table 4. Dust opacity indices and dust disc masses. (1) Source name. (2) Dust opacity index (uncertainty is 0.4). (3) Large grains? 
Y when maximum grain size up to mm-sizes (assumed when /3 < 1), N when /? > 1, ? when 0.7 < /3 < 1.4 due to the uncertainty in j3 
values. (4) Temperature power-law exponent, g, using the IRAS 60 and 100 /im fluxes. (5) Product of dust disc mass and dust opacity 
at 3 mm. (7) The dust opacity evaluated using Eq. |4] (8) The dust disc mass evaluated using the dust opacity Ksmm- (9) The total disc 
mass (gas+dust) assuming a gas-to-dust ratio of 100. 



Sources 


P 


Large (JlRAS M^^^^ x K-^mm K:imm 






P~a-2 


grains? 5 = 2/(3-Qir) IO^'M© x cm'^g"-' cm^g"-' 


10-"' Mq (g/d = 100) 



Chamaelcon 



SY Cha* 






0.3 










CR Cha 


1.3 


7 


0.3 


1.2 


0.4 


3.0 


3.0E-02 


CS Cha 


0.9 


7 


0.4 


1.9 


1.2 


1.6 


1.6E-02 


DI Cha 


1.1 


7 




0.5 


0.7 


0.7 


6.7E-03 


T Cha 


1.1 


7 


0.4 


0.5 


0.7 


0.7 


7.0E-03 


Glass I 


1.2 


7 




0.8 


0.6 


1.3 


1.3E-02 


SZ Cha 


0.9 


7 




0.4 


1.2 


0.9 


9.4E-03 


Sz 32 


1.8 


N 




0.6 


0.1 


4.3 


4.3E-02 


DK Cha 


0.9 


7 


0.7 


12.2 


1.2 


10.6 


l.lE-01 


Lupus 


IK Lup 


0.4 


Y 




0.6 




0.1 


1.4E-03 


Sz 66* 
















HT Lup 


0.0 


Y 


2.9 


2.1 


9.7 


0.2 


2.1E-03 


GQ Lup 


0.2 


Y 


2.1 


0.5 


5.9 


0.1 


8.7E-04 


GW Lup 


0.3 


Y 


0.2 


1.5 


5.6 


0.3 


2.7E-03 


RY Lup 


1.9 


N 


0.7 


0.9 


0.1 


6.9 


6.9E-02 


HK Lup 


0.9 


7 


0.4 


2.3 


1.1 


2.0 


2.0E-02 


Sz 111 


0.5 


Y 


0.2 


1.8 


3.3 


0.5 


5.4E-03 


EX Lup 


0.6 


Y 


0.2 


0.6 


2.3 


0.2 


2.4E-03 


MY Lup 


0.3 


N 


0.4 


1.8 


5.4 


0.3 


3.4E-03 


RXJ1615. 3-3255 


1.6 


N 


0.3 


1.8 


0.2 


7.7 


7.7E-02 



* SY Cha and Sz 66 only have upper limits at 1.2 and 3 mm. 
Boldface sources were resolved at 3 mm. 

Sources with " — " for giRAS have no available IRAS 60 and/or 100 lira fluxes. 



iRicci et al.ll2010al ). Assuming a gas-to-dust ratio of 100, we 
find eight sources have a total disc mass (gas+dust) greater 
th an 0.01 Me, the minim um mass solar nebula according 
to IWeidenschillinS (| 19771 ). and six have a total disc mass 
(gas-|-dust) greater than 0.02 Mq, the minimum mass solar 
nebula according to iHavashil (|l98lf ). 

Note that there is a lot of uncertainty in the dust mass 
calculations. One of the main sources of uncertainty comes 
from the dust opacity, which is dependent o n the chemica l 
composition, size, and shape of the grain (e.g. , lDrainell200^ . 
A second source of uncertainty comes from the assumed dust 
temperature, where a change of ±5 K can cause a ~20% 
change in the dust mass. The uncertainty in the primary 
flux calibrat ion and the uncertainty in the sourc e distances 
of ± 50 pc (jLuhmanl l2007l : iHughes et all 1 1994 IComeronI 
|2008| ) also contributes to the total uncertainty of the dust 
meisses. That said, the method employed here is that gener- 
ally used in the lite r ature (e.g. Andrews fc Williamsl l2005l : 
iLommen et"alll2009l . l2010l : iRicci et al.ll2010al ) and~one can 
expect the disc dust masses to be constrained within a factor 
of 2 to 10 due to these systematic uncertainties. 



4.3 Millimetre-sized grains? 

Our 3 and 7 mm fluxes presented in Table |3] (and in Figs. [6] 
andO indicate that 11 sources do not have a break in their 



spectral slopes at 7 mm, suggesting that thermal dust emis- 
sion is dominant up to 7 mm. From the derived dust opacity 
indices (assuming j3 ~ a-2) in Table U — where the second 
column presents the j3 values and third column indicates 
the existence of grains up to mm sizes — we flnd that 6/11 
sources, all in Lupus, have /3 < 1, suggesting grain growth 
up to mm sizes, while 3 sources have ISM sized grains with 
P > 1. The majority of the Chamaeleon sources (7 sources) 
have 0.7 < /3 < 1.4, making it difficult to analyse the results 
given that the uncertainty in jS is 0.4. 

Assuming the emission from 1-7 mm is optically thin, 
our results would suggest dust grains in Lupus sources are 
generally larger than in Chamaeleon sources, and the lack of 
excess emission observed at 7 mm for Lupus sources could 
suggest less chromospheric emission and potentially more 
evolved systems than Chamaeleon sources and hence po- 
tentially older. However, compared to published ages of the 
different star forming regions, it is unclear if this is indeed 
the case. Lupus 3 and Chamaeleon I have similar ages (3- 
6 Myr), while Lupus 1 is thought to be younger than L upus 
3 (|Hughes et al.lll994l : lLuhmanll2007l : fMerm et al.ll2008l ). In- 
dividual sources wi th signs of large grains have ages rangin g 
from < 1-14 Myrs (|Hughes et al.ill994l : iLawson et al.lll996l ). 
The lack of a break at 7 mm, is an indication of no ex- 
cess emission which could be a signature of age. However, 
no distinction is found in the ages of individual sources 
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Figure 9. Temporal monitoring of 3+6 cm band continuum flux 
at 9 and 5.5 GHz of 3 Chamaeleon sources. The scan lengths 
were 15 minutes with total integration times given in Tabic [BT] 
The fluxes are represented by triangles with corresponding uncer- 
tainties, and the red squares represent 3(t detection values. Note 
that when a source is not detected (fluxes below the red squares), 
UVFIT can sometimes estimates a negative flux. These are included 
for completeness, however the value itself is incorrect. 



with a break (3- 10 Myrs) and sour c es without a break 
(0.9-14.1 Myrs) jHughes et all 1 1994 iLawson et all 1 19961 : 



iHussain et"aLll2009l : ISchisano et al.ll2009h " 

The temporal monitoring results for both mm and cm 
wavelengths presented in Section 13.41 provide further clues, 
allowing us to rule out other sources of emission besides 
thermal dust at 7 mm and beyond. Note of the 7 sources 
with P < 1 only Sz 111 and MY Lup were monitored for flux 
variability at 7 mm. We found no flux variability at 7 mm for 



MY Lup, Sz 111, and conclude that thermal dust emission 
dominates in these sources and hence they have grains up 
to ~1 cm in size. 

Our results also suggest all six sources observed at 
15 mm have excess emission above thermal dust, with tem- 
poral monitoring suggesting the emission comes from ther- 
mal free-free emission. At 3-1-6 cm all three sources were 
found to have some excess emission. Thus it is difficult to 
obtain grain size information for these source beyond 7 mm. 



4.4 Correlating grain growth signatures 

Thus far we have determined signatures of mm grain growth 
in our survey sample by evaluating /3 ~ a - 2, and found 
35% of the discs have grains up to mm sizes. Another 
grain growth signature is the 10 fim silicate feature, which 
probes the warm inner ("-1-5 AU) and upper layers of typ- 
ical T Tauri stars discs, while the mm band probes the the 
cool er outer disc re g ions (> 10 A U) and mid-plane. 

iLommen et~aLl (|2007l . I2OI0I ') suggested that a tentative 
correlation exists between the strength and shape of the 
10 /xm feature and mm spectral slope for a sam ple of discs 
in Cha maeleon, Lupus and Taurus. However, iRicci et al.l 
(|2010bD f ound no such correlati on in their Taurus sample. As 
noted bv lLommen et al] (|2010l ). a correlation between these 
two regions of the disc is unexpected, and if confirmed, it 
could imply that grain growth in the inner upper layers of 
the disc and in the mid-plane occur simultaneously. Here we 
use our sample to further investigate this potential correla- 
tion. 

For this analysis, Spitzer Infrared Spectrograph (IRS j3 
spectrum data from 5.3-16 pm was obtained from the 
Spitzer Heritage Archiv43 for all our sources and an ad- 
ditional seven p Ophiucus sources (to coincide with the 
iRicci et al.ll2010"bl sample), as well as 11 Taurus- Auriga and 
two addi tional Chamaeleon a nd Lupus sources (to coincide 
with the lLommen et al]|2010l sample) - see Table (5) All the 
infrared data presented in Table [5] were processed with the 

S18.18 pipeline 

F ollowing iKessler-Silacci et al. ] (|2005h and lFurlan et al.l 
l|2006h . we define the shape of the 10 ^tm silicate feature 

where 



f 9.8)Jn 



as the fiux ratio at 11.3 /im and 9.8 fj,m, 
Fii.s^im and Fg.g^jm are the integrated fiux over a 0.4 fim 
band centred at 11.3 and 9.8 /xm, and define the strength of 
the 10 ^m silicate feature as: 



strength = 



(F 



10/J.m 



-F, 



t) 



(5) 



where Fio^jm is the integrated fiux from 8-12.4 and Fcont 
is a third-order polynomial fit to the continuum from 5-7.5 
fim and 12.5-16 ^m. 

For each source presented in Table [SI the strength and 
shape of the 10 /im silicate feature were calculated by us- 
ing the three IRS spectrum tables containing the wavelength 



Spitzer Space Telescope is operated by the Jet Propulsion Lab- 
oratory, California Institute of Technolog y under a contract with 
NASA. For more information on IRS see iHouck et al] ||2004|) 

Housed at the Infrared Science Archive (IRSA) at the Infrared 
Processing and Analysis Centre (IPAC), California Institute of 
Technology 
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Figure 10. The 1—3 mm spectral slope a as a function of the 
10 /im strength (top) and shape (bottom) using data from Ta- 
blc[5] Black circles are Cham sources, red ties are Lupus sources, 
blue triangles are p Ophiucus, open diamonds are the isolated 
sources and pink diamonds are Taurus- Auriga sources. The ob- 
tained Spearman rank correlation coefficient r and the corre- 
sponding p-value (probability the obtained r value happen by 
chance) are given for the 10 lira strength and a and the 10 ^m 
shape and a. 



ranges 5-8.6 /im, 9-19 /^m or 13-21 /im, and 7-14 /im (corre- 
sponds to modules SL, SH, LL and SL respectively), with the 
exception of Taurus-Auriga sources which used two tables 
(9-19 /im and 7-14 /im). Note that the method of finding 
the 10 /im strength and shape is sensitive to the wavelength 
range chosen. Here we used the ranges from 5-7.5 /xm, 12.5- 
16 /im, 8-12.4 /im to determine the 10 /im strength, and 
0.4 /im bandwidth to determine t he 10 /im shape, whi ch are 
consistent with the ranges used bv lFurlan et al.l (|2006l ) - see 
Figure [Ell for the third-order polynomial fits to the infrared 
continuum. 

No IRS spectrum was available for SY Cha, DK Cha 
and Sz 111, and only the 9-19 /im table was available for 
GQ Lup and RU Lup and hence these sources are not listed 
in Table [5] T Cha and RXJ1615. 3-3255 are isolated sources 
and were not included in the statistics of the Chamaeleon 
and Lupus sample. SZ Cha has a strong Polycyclic Aro- 
matic Hydrocarbons emission band at 11.3 /im and thus the 
bandwidth ranged of 11.28-11.31 /im was used to determine 

Fll.3,jm. 

The 1-3 mm spectral slope s a are from Tab le [3] 
for Chamaeleon and Lupus, from iRicci et al. (l2010b|) fo r 
the p Ophiucus sour c es and from iRodmann et al.l ~ 20061 ) ; 
[Andrews fc WiUiamd (|2007bl) for Taurus-Auriga - see Ta- 
ble (H 



Table 5. Strength and shape of the 10 /tm silicate feature for all 
Chamaeleon and Lupus sources from this work, and supplemen- 
tary p Ophiucus and Taurus-Auriga sources. (1) Source name. 
(2) 10 /im strength values. (3) 10 /im shape. (4) 1—3 mm spectral 
slope. 



Source 


Strength 


Shape 


a 


Chamaeleon 


CR Cha 


1.07 


0.71 


3.4 


CS Cha 


0.40 


0.58 


2.9 


DI Cha 


0.42 


0.84 


3.2 


Glass I 


1.19 


0.72 


3.1 


SZ Cha* 


0.60 


0.80 


2.9 


Sz 32 


0.18 


1.05 


3.8 


WW Cha" 


0.77 


0.78 


2.8 


Lupus 


IK Lup 


0.13 


0.83 


2.4 


HT Lup 


0.31 


0.85 


2.0 


GW Lup 


0.28 


0.92 


2.3 


RY Lup 


1.24 


0.60 


3.9 


HK Lup 


0.67 


0.73 


2.9 


EX Lup 


0.54 


0.74 


2.5 


MY Lup 


0.55 


0.81 


2.3 


IM Lup" 


0.51 


0.85 


3.6 


Isolated 


T Cha 


-0.03 


0.77 


3.1 


RXJ1615. 3-3255 0.93 


1.02 


3.6 


p Ophiucus 


SR 4 


0.33 


1.38 


2.5** 


IRS 49 


0.51 


0.98 


1.8^* 


WSB 52 


0.06 


1.01 


1.8* 


WSB 60 


0.28 


0.99 


1.9'' 


DoAr 44 


1.39 


0.65 


2.2'' 


RNO 90 


0.39 


0.79 


2.3'' 


Wa Oph 6 


0.18 


0.90 


2.4^* 


Taurus-Auriga 


DC Tan 


-0.03 


0.67 


2.1= 


DO Tau 


0.18 


1.24 


2.3= 


AA Tau 


0.32 


0.54 


3.2^* 


CI Tau 


0.54 


0.58 


2.2'' 


DL Tau 


0.10 


0.61 


2.0'' 


DM Tau 


0.41 


0.80 


2.9'' 


DN Tau 


0.14 


0.62 


2.3'' 


DR Tau 


0.22 


0.54 


2.2'' 


FT Tau 


0.20 


0.61 


1.8'' 


CM Aur 


0.84 


0.48 


3.2'' 


AS 205 


0.50 


0.54 


0.7'' 



* SZ Cha has a strong Polycyclic Aromatic Hydrocarbons 
emission band at 11.3 /im. 

" Additional Cha maeleon and Lupus sources from 
iLommen et al.l 1 I2OIOI). 



a obtained from 
= a obtained from 
a obtained from 



Ricci et alj bOlObI). 
Rodmann et al.l i200p ). 



Andrews fc Williamsl l|2007bh 
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The millimetre spectral slope a as function of 10 
strength and shape where plotted in Fig. 1101 To test 
these two correlations we used the statistical package R 
||R Development Core Teamll201ll ) to determine the Spear- 
man rank correlation coefficient, where the rank coefficient, 
r, is a measure of the how two values are related (r = weak 
correlation, r = ±1 strong correlation), p- value is a measure 
of the probability the obtained r value happen by chance. 
With a p- value of 0.17 and 0.39, we found the rank coef- 
ficients for the 15 sources in Chameleon and Lupus to be 
r = 0.37 and r = -0.24 for strength versus a and shape ver- 
sus a respectively - see Table [H] Thus we conclude that 
the correlations between 10 ^m strength versus a and the 
10 ^m shape versus a cannot be confirmed with this sam- 
ple. Spearman statistics was then performed on the sample 
of 36 sources from the four star forming regions and on each 
individual star forming region listed in Table [5] 

The results suggest no correlation is present for each in- 
dividual star forming regions (Lupus 1 and 3 only have three 
sources, thus were not included in the individual analysis) - 
see Table El A weak correlation between the strength of the 
10 fim silicate feature and a is found for the whole sample, 
while no correlation appears present between the shape of 
the 10 fim silicate feature and a. 

There are multiple processes affecting the shape of the 
10 fim silicate feature (grain growth, crystallisation and 
opacity) that would make the detection and interpreta- 
tion of a possible correlation difficult. For example both 
RXJ1615. 3-3255 and WSB 52 have a 10 fim shape of 1.02 
and 1.01 respectively, suggesting a flat shape structure for 
both sources. However, this is not the case RXJ1615.3- 
3255 has a strong, boxy 10 fim feature, while WSB 52 
has a very fiat structure - see Fig. lEll It was only when 
Spitzer allowed for sufficient statistics that such ambiguities 
were cancelled out and the correlation between the strength 
and the shape of the 10 fim feature was confirmed (e.g. 
iKessler-Silacci et al.ll200^). 



Kessler-Silacci et al 



(|2006l ) have shown a similar corre- 
lation exists between the 10 and 20 /im strengths, with some 
outliers with strong 10 fim features but weak 20 nm features. 
To reproduce these results, they had to increase the size of 
the grains emitting at 20 fim to a minimum radius of 2-3 fim, 
implying the emission at 20 ^m is from a region in the disc 
deeper than the grains emitting at 10 nm with a minimum 
radi us of 0.1 fim. A r ecent study performed on brown dwarfs 
by (|Riaz et al.ll20lj ) has found the flattening and shift of 
the peak position of the 10 fim silicate feature to have no 
correlation with an increase in the grain sizes or a higher 
degree of crystallinity in the disc. The four outliers identi- 
fled in the study have a mix of small, large and crystalline 
grains, with a weak 20 fim. Other studies have found similar 
outliers in T Tauri stars which were ex plained by enhanced 
dust s ettling in the out er disc region jKessler-Silacci et al.l 
l2007l ). lRiaz et all (|2012l ) suggest that the presence of these 
outliers indicates that the timescales for dust processing are 
different for the inner and outer regions of discs. 

Compared to previous analysis of this correlation, our 
sample is the largest (36 sources co mpared to 7 sou r ces in 
iLommen et al.l (12007 ), 30 sources in iLommen et al.l (|2010l ) 
and 21 sources in iRicci et al.l (|2010bl ')'). excludes sources for 
which only upper limits on a are known, contains no upper 
limits for 1-3 mm spectral slope a, and all the IR data were 



process and analysed consistently. However, the sources used 
may have several intrinsic differences. In particular, Cha I 
and p Ophiucus are very different e nvir onments, with Cha I 
sources being fairly isolated (e.g. iLuhman ,200^) while p 
Ophi ucus sources are deeply embedded (e.g. iMcClure et al.l 
I2OI0I ). Although the p Ophiucus values in Table [S] were 
no t corrected f or exti nction, similar results were obtained 
bv lRicci et al] (|2010bl ) using extinction corrected values for 
the 10 fim shape. Ideally we would like the same sample 
number per star forming region for this type of analysis. 
A larger sample may also mitigate the intrinsic differences 
in the 10 ^m silicate feature, allowing for a more deflnite 
conclusion on the presence or absence of a correlation with 
a. 



5 CONCLUSIONS 

Continuum observations were carried out with ATCA at 
3, 7, 15 mm, and 3+6 cm for 20 T Tauri stars located in 
Chamaeleon and Lupus. We analysed the mm fluxes in or- 
der to determine the spectral slopes, maximum grain sizes, 
and dust masses in these discs. Using supplementary data 
from the literature we conducted temporal monitoring of the 
fluxes of a subsample of our sources over short (less than 
one day) and long (months to years) timescales to help con- 
strain the emission mechanisms present in these discs. We 
also analysed the potential correlation between the millime- 
tre spectral slope and the st rength and shape of the 10 fim 
silicate feature suggested bv ILommen et all (I2OO7I . I2OIOI '). 

(i) Our 3 and 7 mm continuum fluxes show that 11 
sources do not have a break in the spectral slope at 7 mm, 
suggesting thermal dust emission is dominant to wave- 
lengths as long as 7 mm. We found that 6 of those sources 
(all in Lupus) have a dust opacity index less than unity (as- 
suming /? ~ Q - 2), suggesting grain growth up to at least 
mm sizes. 

(ii) We obtained dust disc mass estimates ranging from 
10^^-10"^ Me. Assuming a gas-to-dust ratio of 100, we flnd 
eight sources have a total disc mass (gas+dust) greater 
th an 0.01 Me, the minim um mass solar nebula according 
to IWeidenschillind (| 19771 ): and six have a total disc mass 
greater th an 0.02 Me, th e minimum mass solar nebula ac- 
cording to iHavashil (|l98lh . 

(iii) All six Chamaeleon sources observed at 15 mm, have 
excess emission above thermal dust. At 3-f6 cm, DK Cha, 
T Cha and Sz 32 have flux upper limits which suggest the 
emission at 3+6 cm is due to an ionised wind and/or chro- 
mospheric emission. 

(iv) The monthly and yearly temporal flux monitoring 
revealed no excess emission at 7 mm for MY Lup, Sz 111 
and CS Cha, indicating thermal dust emission dominates in 
these sources and hence that they have grains up to ~1 cm 
in size. At 15 mm the excess emission for CR Cha, CS Cha, 
DI Cha, T Cha, Sz 32 and DK Cha is not from a fast varying 
chromospheric emission, but most likely from thermal free- 
free emission consistent with the spectral slope of 0.6. 

(v) Supplementing our sample with 16 other sources from 
Taurus-Auriga and p Ophiucus, no correlation between the 
shape of the 10 ^m silicate feature and the 1-3 mm spectral 
slope a was found, while the strength of the 10 pm feature 
appears to correlate weakly with a. 
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Table 6. Spearman rank correlation coefficient values and percent confidence levels (p- values) for each star forming region in Table [5] 
used in the correlation analysis for the 10 silicate feature strength and shape with the 1-3 mm spectral slope a. 







Strength 


vs. a 


Shape 


vs. a 




# of sources 


rank coefT. 


p- value 


rank coeff. 


p- value 


Cham & Lupus 


15 


0.37 


0.17 


-0.24 


0.39 


All 


36 


0.40 


0.02 


0.02 


0.90 


Cha I 


7 


-0.22 


0.64 


0.36 


0.43 


p Ophiucus 


7 


0.02 


0.97 


-0.04 


0.93 


Taurus 


11 


0.32 


0.34 


-0.05 


0.87 



ALMA will be vital in obtaining high resolution maps 
at high sensitivities of protoplanetary discs in the southern 
hemisphere at (sub)mm wavelengths, providing a better un- 
derstanding of disc sizes and further analysis of mm grain 
growth signatures. However, at 7 and 15 mm, the ATCA 
and VLA will continue to provide invaluable information on 
other emission mechanisms present in protoplanetary discs. 
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APPENDIX A: COMMENTS ON INDIVIDUAL 
SOURCES 

Al Chamaeleon 

Chamaeleon comprises of three clouds - C ha I, II and HI 
- located at a distance of ~150-160 pc jLuhmanl l2008l : 
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iFurlan et allbOOgl '). This survey focused on T Tauri stars on 
the Cha I cloud (with the exception to DK Cha in Cha II). 
Cha I has two subclusters, with ages ranging from 3-4 and 
5-6 M yr in the sout hern and northern subcluster respec- 
tively (|Luhmanll2007l '). 

CS Cha is a cl ose binary with a < 0.1" separation 
l|Espaillat et all 120071 '). A near-infrared flux deficit was de- 
tected, suggestive of an inner disc hole, modelled to be 
~43 AU having ~ 5 x 10~^ lu n ar masses of dust from 0.1 
to 1 AU (|Espaillat et al.ll2007h . iKim et al.l (|2009l ) observed 
an excess at 5-8 ^.m higher than typical T Tauri stars. 

DI Cha is a binary system w ith ~ 4". 6 separation 
l|Ghez et al.ll 19971 : ICorreia et "alll2006l ). The binary is yet to 
be resolved in (sub) mm wavelengths (|Henning et al.|[l993l : 
iLommen et al.ll2007l ). 

T Cha is an isolated source. iGeers et all \2mi ) de- 
tected a weak 11.2 /im feature and a marginal 3.3 /xm fea- 
ture. Both features were unresolved, co rresponding t o a dis - 
tance > 13 and < 22 AU respectively. iBrown et all (|2007l ') 
detected a possible gap in the disc at 66 AU. 

Glass J is a binary with separation ~ 2 ". 5 ± 
0.5 llFeigelson fc Krisd Il989l: iReipurth fc Zinnecked 1 19931 : 



CheUi et all 



iGhez et al.lll997l: lLuhmanll2004l ). Glass la is a K4 WTTS 



198S) and Glass lb an IR G-type source 
Feigelson fc Krisd 1 19891 '). The 10 /im feature was as- 
signed to Glass lb. no detection was obtained for la 
ijStanke fc ZinneckejboOOl ). The binary system has not been 

resolved at mm wavelengths. 

SZ Cha was classified by iManoi et al.l l|201ll ) as an 
object with an inner and outer optically thick disc sepa- 
rated by optically thin gap with the 10 lira silicate feature in 
emission. It is an accreting source (|Luhmanll2003 ) with two 
potential companions 2M J05 81804-777197 at 5 ". 2 and 2M 
J10581413-7717088 at 12". 5 (|Ghez et al.lll997h . which are 
not confirm ed members of the Cha I associatio n jLuhmanl 
I2OO7I . l2008h . The IRS spectra from iKim et all (|2009l l de- 
tected polycyclic aromatic hydrocarbon emission at 11.3 /im. 

DK Cha is the br i ghtes t IR and submm source 
in Cha II (jWhittet et al.l 1997l l. found to be transition- 
ing from an emb edded disc to a proto planetary disc 
through modelling (jvan Kempen et al.ll2O10l ). and still sur- 
rounded by a protostellar enve lope and molecular outflow 
l|van Kempen et aLllioO^ . |2009| ) . 



A2 Lupus 

The Lupus southern star for ming region is co mposed of 8 
clouds (Lupus 1 trough 8) (jHara et al.l 1 19991 ). This work 
focused on sources in Lupus 1, 3 and 4 at dis tances of 
150±20 pc, 200±20 pc and 165±15 pc respectively (|Comer6nl 
l2008h . 

IK Lup and S z 66 axe a binary pair with a 6". 4 
(960 AU) separation jSciiwartjl 19771 : iReipurth fc Zinneckeil 
ll993l : lMonin et aLlliooel ). Sz 66 was not detected f or the 1.2 
and 3.2 mm observations bv lLommen et all ()2010h . in both 
cases the beam size was sufficient to resolved the binary 
system. Sz 66 was de tected with Spi t zer an d has a 10 fj,m 
silicate feature - see iLommeri et all (|20ld ). However, the 
lack of cold dust from 1.2 to 7 mm would suggest there is 
no circumbinary disc, and if a cold dust disc is present its 
mass is too low to have been detected by our observations. 

HT Lup is a triple system (Chez et al., 1997), with 



one companion at 2". 8 ± 0.1 (Reipurth fc Zinnecker 1993; 
Brandner et al., 1996) and another at 0".107±0.007 (Ghez et 
al. 1997). Heyer & Graham (1989) found a reflection nebula 
near HT Lup and Ha emission possibly associated with an 
HH jet. 

GQ Lup has a known companion, GQ Lup B. at a dis- 
tance of 0.7", which may be a brown dwarf JCuenther et al 



20051: i Neuhauser et al]|2008ar ) or a planet (jNeuhauser et al 



2008bl ). The system was observ e d at 1 .3 mm with the Sub- 



millimetre Array by iDai et al] (|201[i ) who did not detect 
thermal dust emission from QQ Lup B and model a com- 
pact disc of mass ~Mjup around GQ Lup, which is consistent 

with our disc mass estimate. 

Sz 111 was identified as a cold disc by iMerm et al.l 

(j2008h who suggests the fluxes from the IRAC and MIPS 
band 1 were more representative of a debris disc with no 
gas present (where OK-24fim < -1-6), while showing an ex- 
cess at 70 /xm t ypical of a classical T T auri star. Cold discs 
were defined bv lCalvet et all |200i) and lBrown et al.l (|2007l ) 
as optically thick discs with large inner holes, where poten- 
tial large grain growth ma y occur. Sz 111 i nner hole was 
determined to be ~ 70 AU jMen'n et al.ll2008h . 

EX L up is a highly variable M-type star with no known 
companion (jHerbidBoOTl : iKospal et al.ll2008l l. It is located 
at th e edge of a gap between Lupus 3 and 4 (jCambresvl 
I1999I). It has a quiescen t accretion rate of ^ 4 x lO'^^Mo/yr 
(jNatta fc Testili2004 : Sipos et al.ll2009h . thou gh five flares- 
ups have been observed from 1995- 2 008 dHerbig et al 



2OOII: iHerbid l2007l : iGrosso et all |2010|: lAspin et al l 201 
Goto et al.ll201lf l. In February of 2008. lKospal et all (| 20oi 



observed the extreme outburst of EX Lup. They found there 
a wealth of metallic lines with no added absorption features, 
which lead to the conclusion the brightness increase was due 
to an inc rease in accretion and possibly stellar winds. Mod- 
elling bv lSipos et al.l (|2009l ) found the 10 /xm feature could 
be recreated with amorphous silicates with olive and pyrox- 
ene stoichiometry, and the disc an inner and outer radius 
of 0.2, 150 AU respectively. They concluded the inner hole 
may be the cause of the outburst, however more evidence is 
required to confirm this hypothesis. 



APPENDIX B: OBSERVING DETAILS 

The Australia Telescope Compact Array (ATCA) was used 
to observed 9 sources in Chamaeleon and 11 source in Lupus 
star forming regions at 3, 7, 15 mm and 34-6 cm. A log of 
the ATCA observations is given in Table [Bl] 



APPENDIX C: RESULT DETAILS 



The complete results used in Section 13.11 analysis to deter- 
mine if a source was extended are given in Table ICll and 
the full spectral energy distribution plots are given in Fig- 
ures |Cl] and jCl 
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Table Bl. ATCA observation log. (1) Observation date. (2) Sources observed. (3) Total integration time used in analysis in this paper. 


(4) Central frequency pair of observations. (5) Array configuration*. (6) Notes. 


Observation dates Sources 


Tint 


Frequency pairs Array Notes 






(min) 


(GHz) Config. 


May-OQ" 


SZ Cha 


200 


43, 45 H214 




CS Cha 


80 






CR Cha 


80 






Class T 


150 






Sz 32 


130 






DI Cha 


100 






T Cha 


60 




May-Og' 


Sz 65-Sz 66 


70 


43, 45 H214 




HT Lup 


40 






GW Lup 


60 






RY Lup 


60 






EX Lup 


60 






MY Lup 


50 






Sz 111 


60 






RXJ1615.3-3255 


50 






HK Lup 


50 




May-lO'" 


Glass I 


50 


93, 95 H214 




DI Cha 


104 






Sz 32 


180 




Jul-10" 


DI Cha 


80 


EW352 




DK Cha 


43 






CS Cha 


60 




Jul-lO" 


DK Cha 


54 


43, 45 EW352 


Aug-lO' 


EX Lup 


70 


93, 95 H168. Antenna 2 oflSine. 




GO T,ur> 


44 






ATCA160043-415442 


73 




Aug-10* 


RXT161,5 3-32,55 


60 


43, 45 H168 




GQ Lup 


80 






SY Cha 


30 




Jul- 11= 


DK Cha 


18 


17, 19 H214 Errors with antenna 6 polarization. 




Sz 32 and WW Cha 


42 






T Cha 


39 






DI Cha 


100 






CR Cha 


42 






CS Cha 


105 




Jul- 11"= 


DK Cha 


180 


5.5, 9 EW352 July 28 - antenna 5 oflJine. 








July 29 - antenna 3 offline for last hour. 








Flux calibrator bootstrapped from gain calibrator. 




Sz 32 and WW Cha 


125 






T Cha 


170 












* ATCA array 


configuration information: |http:/ /www.narrabri.atnf.csiro.au/operations/array_configurations/upcoming_configs.html| 


" Phase calibrator QSO B1057-797, primary flux 


calibrator Uranus. 



Phase calibrator QSO B1600-44, primary flux calibrator Uranus. 

Phase calibrator QSO B1057-797, primary flux calibrator QSO B1934-638. 



to 

Table CI. Summary of 3 and 7 mm flux fittings. (1) Source name. (2), (3) Point source and Gaussian fit fluxes. (4) Gaussian size obtained from the Gaussian fit. (5) Synthesiscd tf^ 

beam size using natural weighting. (6) Image RMS (natural weighting). (7) Factor of a the point source fit is below the Gaussian fit (Fg=Fp+n(T). (8) Description of the amplitude as 

a function of u-v distance plots (UVAMP) presented in Figs. [2] and |3] (9) If the emission was considered resolved or not. 



Sources 


Fp 


Fg 


Gaussian size 


Beam size 


RMS 


n 


UVAMP" 


Resolved'^ 




fmjv) 


f mjv) 


(arcsec) 


(arcsec) 


(mJy /beam) 


Fg=Fp-|-nCT 


F,D,S 


Y,N 








Q 
O 












CS Cha 


7.8 ± 0.8 


9.4 ± 0.9 


14.3 ± 1.9 X 1.8 ± 0.2 


13.2x1.4 


0.6 


2.7 


D 


Y 


DI Cha 


2.3 ±0.5 


2.3 ±0.4 


4.2±1.4x 1.9 ±0.6 


4.6x1.8 


0.3 


0.0 


F 


N 


Glass I 


4.1 ±0.1 


4.4 ±0.1 


2.5±0.2 X 2.0 ±0.1 


2.4x1.8 


0.1 


3.0 


F 


N 


Sz 32 


3.1 ±0.3 


3.1 ±0.2 


2.4±0.6 X 1.9 ±0.5 


2.4x1.8 


0.2 


0.0 


F 


N 


DK Cha 


33.0 ±3.0 


49.8 ±1.3 


7.7 ±0.5 X 2.4 ±0.2 


7.1x1.3 


1.3 


12.9 


D 


Y 


GQ Lup 


3.6 ±0.3 


3.5 ±0.8 


2.9±0.9 X 1.4 ± 1.1 


2.8x2.0 


0.4 


-0.3 


F 


N 


EX Lup 


2.0 ±0.2 


1.9 ±0.3 


3.2±0.4x 1.7±0.2 


2.9 x 1.9 


0.3 


-0.3 


F 


N 








7 


mm 










Chamaeleon 


SY Cha 


<1.0= 






16.0x5.4 


0.3 








CR Cha 


1.7 ±0.1 


1.8±0.1 


5.1±0.4x4.8±0.4 


5.0x4.3 


0.1 


1.0 


D 


N 


CS Cha 


1.5 ±0.1 


1.7±0.1 


6.0 ±0.6x4.8 ±0.5 


5.0x4.3 


0.1 


2.0 


D 


N 


DI Cha 


0.9 ±0.1 


0.9±0.1 


4.7±0.5 X 4.4 ±0.4 


4.8x4.2 


0.1 


0.0 


F 


N 


T Cha 


3.0 ±0.1 


3.0±0.1 


5.2±0.2 X 4.1 ±0.1 


5.1x4.1 


0.1 


0.0 


F 


N 


Glass I 


0.7 ±0.1 


0.7±0.2 


5.2±0.9x4.5±0.8 


4.9x4.1 


0.1 


0.0 


F 


N 


SZ Cha 


0.7 ±0.1 


0.8 ±0.1 


5.8±0.9x4.7±0.8 


5.0x4.1 


0.1 


1.0 


F 


N 


Sz 32 


1.2 ±0.1 


1.2±0.1 


5.3±2.0x3.8±1.4 


5.0x4.0 


0.2 


0.0 


S 


N 


DK Cha 


6.6 ±0.6 


8.4±0.2 


54.9 ±4.5x3.1 ±0.3 


47.4x2.5 


0.7 


2.6 


D 


Y 


Lupus 


Sz 65-66 


0.9±0.7 


1.0 ±0.1 


4.9±0.5x3.5±0.4 


4.5x3.5 


0.1 


1.0 


F 


N 


HT Lup 


3.4 ±0.1 


3.5±0.1 


4.6±0.1x3.5±0.1 


4.4x3.5 


0.1 


1.0 


F 


N 


GQ Lup 


0.6 ±0.1 


0.6 ±0.1 


4.2±1.3x 1.1 ±2.8 


6.7x4.4 


0.1 


0.0 


F 


N 


GW Lup 


0.6 ±0.1 


0.6 ±0.1 


5.2±1.0x3.6±0.7 


4.8x3.5 


0.1 


0.0 


F 


N 


RY Lup 


1.0 ±0.1 


1.0 ±0.3 


4.9±0.3x3.5±0.2 


4.7x3.5 


0.1 


0.0 


F 


N 


HK Lup 


1.0 ±0.1 


1.0±0.1 


5.1±0.5 X 3.4 ±0.3 


5.2x3.5 


0.1 


0.0 


F 


N 


Sz 111 


0.5 ±0.1 


0.5 ±0.1 


5.8±1.4x3.3±0.8 


5.2x3.5 


0.1 


0.0 


F 


N 


EX Lup 


1.4 ±0.1 


1.4±0.1 


5.1±0.5x3.5±0.9 


4.8x3.5 


0.1 


0.0 


F 


N 


MY Lup 


1.1 ±0.1 


1.1 ±0.3 


4.8±0.7x3.7±0.5 


4.8x3.4 


0.2 


0.0 


S 


N 


RXJ 1615.3-3255 


0.8 ±0.2 


0.8 ±0.2 


7.2 ±4.4x6.2 ±4.9 


9.5x4.1 


0.2 


0.0 


D 


N 



" UVAMP: F - flat, D drop, S - sinusoidal. 
* Resolved: Y - yes, N — no. 
A 3(T upper limit is given for non-dctcctions. 
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APPENDIX D: MONITORING DETAILS 

The complete results used for the monitoring analysis pre- 
sented in Section r3.4l The 7 mm fluxes are given in TablelDll 
f5 mm in Table |D2] and 3+6 cm in Table lD3l 



APPENDIX E: STRENGTH AND SHAPE 

The third-order polynomial fits to the infrared continuum 
obtained from the Heritage Infrared Archive of each source 
are present ed. Spitzer data has been published in previ- 




iLommen et al.ll2010l l. 
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Figure CI. Full S EDs for Chamaeleon sources. Data obtained from the literature and this work. The stellar photosphere was obtained 
using the method of lFurlan et al.l l|2006l ) normalising sources of spectral type later than G by their J-band flux, and earlier spectral types 
by their V-band flux. The photometry data were obtained from published catalogs. Legend: * WFI (Wide Field Imager); + 2MASS (Two 
Micron All Sky Survey); x CTIO (Cerro Tololo Inter-American Observatory); v USNO (U.S. Navel Observatory); ■ IRAC (Infrared 
Array Camera on Spitzer Space Telescope); o IRAS (Infrared Astronomy Satellite); o MIPS (Multiband Imaging Photometer for Spitzer); 
• ESO-TIMMI2 (Thermal Infrared Multi Mode Instrument 2) o n the ESO 3.6-m telescope at La Silla; A LABOCA (LA rge BOlometer 
CAmera) on the 12-m APEX telescope jBelloche et al.l 2011hl ); A SEST (Swedish-ESO 15m Submillimeter Tele scope). iHenning et alj 
lll993l 'l: □ SMA (Submillimeter Array'). iLommen et al.l l 201Cll 'l: ATCA (Australia Telescope Compact Array). iLommen et al.l l|2007l . 
|201C|) and this work. 
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Figure C2. Same as Fig. ICll but for Lupus sources. Note the stellar photosphere of MY Lup was normalised to the J-band flux since 
its spectral type is unknown. 
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Table Dl. Results of 7 mm flux monitoring of RY Lup, Sz 111, MY Lup, RXJ1615.3-3255, CS Cha and Sz 32 from this work and 
the Hterature. (1) Source name and dust opacity index /3 ( see Section 14. l|l . (2) D ate o f observation in YYY YMMDD format. (3) ATCA 
array configuration". (4) References: pre-CABB data from lLommen et alj | |2009|) and iLommen et al. 1 1 I2OIOI) . and CABB data from this 
work. (5), (6) continuum fluxes from point source flts obtained at approximately 6.7 and 7.0 mm respectively with the RMS given in 
parenthesis in units of mjy/bcam. 



Source 


Date 


Array 


Reference 


F(6.7 mm)*" 


F(7 mm)*" 










(mJy) 


(mJy) 


RY Lup 


20080801 


H214 


Lommen et al. (2010) 


— 


< 0.6(0.2) 


/3 = 1.9 


20090529 


H214 


This work 


< 1.2(0.4) 


< 0.9(0.3) 




zuuyuoou 




This work 


1 ^ 4- c\ ^ (c\ ^\ 

±, i ± U. i I^U. ± ) 


U.O ± u. 1 l^U. i ) 


MY Lup 


20080805 


H214 


Lommen et al. (2010) 


1.3(0.1) 


— 


P = 0.3 


20090529 


H214 


This work 


< 1.2(0.4) 


< 1.2(0.4) 




^uuyutjou 


H214 


_L lliD WUliS.. 


1 n -1- n 1 ("n 1 1 




RXJ1615.3-3255 


20080801 


H214 


Lommen et al. (2010) 


< 0.5(0.2) 




/3 = 1.6 


20100821 


H168 


This work 


0.8 ±0.2(0.2) 


0.8 ±0.2(0.2) 


Sz 111 


20080801 


H214 


Lommen et al. (2010) 


< 0.6(0.2) 




P = 0.5 


20090529 


H214 


This work 


< 0.9(0.3) 


< 0.9(0.3) 




20090530 


H214 


This work 


0.3 ±0.1(0.1) 


0.3 ±0.1(0.1) 


CS Cha 


20080426 


6A 


Lommen et al. (2009) 




1.0±0.3(0.1) 


13 = 0.9 


20080705 


1.5B 


Lommen et al. (2009) 


< 0.8(0.3) 


< 0.7(0.2) 




20080706 


1.5B 


Lommen et al. (2009) 


< 1.1(0.4) 


1.4±0.3(0.2) 




20090528 


H214 


This work 


1.7±0. 4(0.3) 


1.2±0.3(0.3) 




20090529 


H214 


This work 


2.0±0.3(0.3) 


1.9±0.2(0.2) 


Sz 32 


20080331= 


H168 


Lommen et al. (2009) 




0.8±0.1(0.2) 


13 = 1.8 


20090528 


H214 


This work 


< 0.5(0.2) 


< 0.5(0.16) 




20090529 


H214 


This work 


1.5±0.5(0.3) 


1.0±0.4(0.2) 



" ATCA array conflguration: |http: / / www.narrabri.atnf.csiro.au/operations / array-Configurations / upcoming_configs.html| 
A 3(T upper limit is given for non-detections. 

Sz 32 was detected in the fleld of view of WW Cha observations. 
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Figure El. The third-order polynomial fit to the infrared continuum used to obtain the strength of the lOfim silicate feature. The 
Spitzer infrared data were obtained from the Heritage Archive. 
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Figure El. Continued. 
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Table D2. Results of the 15 mm flux monitoring of 6 Chamaeleon 
sources, listing the highest and lowest fluxes obtained from a point 
source flt in the u-v plane. A 3cr value is given for non-detections. 
RMS = 0.1 mjy/beam. 







17 GHz 


19 GHz 


Sources 


Time 


Fhigh 


-^low 


Fhigh 


Flow 




(minutes) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


CR Cha 


3 


0.6 ±0.1 


< 0.3 


0.6 ±0.2 


< 0.2 


CS Cha 


5 


0.5±0.1 


< 0.3 


0.4 ±0.1 


< 0.2 


DI Cha 


5 


< 0.3 


< 0.3 


0.3±0.1 


< 0.3 


T Cha 


3 


0.9 ±0.2 


< 0.3 


0.5 ±0.2 


< 0.3 


Sz 32 


5 


0.8±0.1 


< 0.2 


0.7±0.1 


< 0.2 


DK Cha 


3 


2.0 ±0.2 


1.5 ±0.2 


1.8 ±0.2 


1.3 ±0.2 


Table D3. 


Results of 


the 3+6 


cm flux 


monitoring 


of 3 



Chamaeleon sources, listing the highest and lowest fluxes ob- 
tained from a point source fit in the u-v plane. A 3cr value 
is given for non-detections. RMSs.sqhz =0.1 mJy/beam and 
RMSgGHz £ 0-2 mJy/beam. All sources had a scan length of 
15 minutes. 





5.5 GHz 


9 GH 


z 


Sources 


Fhigh 


-^low 


Fhigh 


Flow 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


T Cha 


0.9 ±0.1 


< 0.3 


0.5±0.1 


< 0.5 


Sz 32 


0.8 ±0.1 


< 0.3 


0.7 ±0.1 


< 0.5 


DK Cha 


0.7±0.1 


< 0.3 


1.1 ±0.1 


< 0.4 



